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ASYMMETRIC DIAMOND ISING-HUBBARD 
CHAIN WITH ATTRACTION 



The ground state and thermodynamic properties of an asymmetric diamond Ising-Hubbard 
chain with the on-site electron- electron attraction has been considered. The problem can be 
solved exactly using the decoration-iteration transformation. In the case of the antiferromag- 
netic Ising interaction, the influence of this attraction on the ground state and the temperature 
dependences of the magnetization, magnetic susceptibility, and specific heat has been studied. 
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1. Introduction 

The spin-chain problem, which is solved exactly with 
the use of the decoration-iteration transformation [1- 
y| attracts interest, because it allows certain features 
in the properties of complicated spin systems and 
magnetic materials to be studied [5|-l7| . These are the 
intermediate plateaux on magnetization curves and 
extra maxima on the temperature dependence of the 
heat capacity. The spin chains also enable the interre- 
lation between a geometrical frustration and quantum 
fluctuations to be analyzed |8l4ld|. Moreover, they 
can serve as models for the quantitative description 
of magnetic properties of materials [7|,lll|. That is 
why the one-dimensional models, which can be solved 
exactly with the use of the decoration-iteration trans- 
formation, are actively studied |l2l42l| . 

Work [16], in which the properties of an asym- 
metric diamond Ising-Hubbard chain were considered 
without retard for the on-site electron-electron inter- 
action, has started the researches of exactly solved 
(by applying the decoration-iteration transformation) 
Ising-Hubbard systems [22|, |23|. This chain reveals 
the following features: the and 1/3 magnetization 
plateaux [ig, one [16| or two [2J| additional low- 
temperature peaks in the zero-field heat capacity, and 
the considerable adiabatic magnetocaloric coefficient 



17| . The influence of the on-site Coulomb repulsion 
of electrons on the property of this chain was stud- 
ied in work 2J|. In particular, in the interval with 



substantial repulsion, the zero-field heat capacity was 
shown to have an additional maximum at high tem- 
peratures. 

If the chain of electrons interacts with local 
phonons, an effective on-site electron-electron attrac- 
tion can be realized in it 25l427|. By analogy, in this 
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work, the properties of an asymmetric diamond Ising- 
Hubbard chain [16] with the on-site electron-electron 
attraction is considered. In particular, in the case of 
the Ising antiferromagnetic interaction, when a geo- 
metrical frustration in the chain emerges, the influ- 
ence of this attraction on the ground state and ther- 
modynamic properties will be analyzed. 

2. Model and Its Exact Solution 

Consider an asymmetric diamond Ising-Hubbard 
chain |l6l.l24j with the on-site electron-electron at- 
traction (Fig. 1) in a magnetic field. The primitive 
cell of this chain is determined by the k-th and (fc-|-l)- 
th nodes, both occupied by the so-called Ising spins, 
/ife, which are coupled with neighbor spins by means 
of the Ising interaction. Two electrons with the on- 
site attraction between them execute quantum jumps 
over two interstitial positions (fc, 1) and (fc, 2) in the 
primitive cell. 
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Fig. 1. Schematic diagram of a fragment in an asymmetric 
diamond Ising— Hubbard chain. SoUd and hollow circles de- 
note nodes and interstitial positions, respectively. For the fc-th 
primitive cell, the nodal spins, /ij,, the 2-components of total 
spins in the interstitial positions, Sk,i, and the parameters of 
the Ising interaction for bonds along the diamond sides, li, are 
indicated 

The Hamiltoniaii % of the chain is written down as 
a sum of ceh Hamihonians Hk, 



N 



■H=^nk, 



k=l 



2 



The replacement of the on-site electron-electron 
attraction by a repulsion in Hamiltonian JT]) trans- 
forms it into the chain Hamiltonian from work [2J|, 
which was solved exactly by applying the decoration- 
iteration transformation. Therefore, the exact solu- 
tion for our chain can be obtained by replacing the 
on-site electron-electron repulsion by the attraction 
{U -^ -U) in aU results of work [24]. The spec- 
trum of the Hamiltonian Hk obtained in such a way 
looks like 

EiilJLk, tJ-k+i) = — 7^ — {^J■k + fJ.k+i)-hc- 



-y(/"fe+M/c+i), 



£2{p-k,lJ-k+l) = 7:, (M/c +M/c+l) + /ie 



-y(Aifc+Mfe+i), 



hi. 



£3{^Ik,^^k+l) = Ail^fc - Mfe+il - -^ifj-k + ^j-k+i), 



■ fJ-killSk.l + l2Sk,2) + Mfc+l(^2<5'fe,l + hSk, 



hciSk,i + Sk^2) ~ -7:h\{iik + A^fe+i), 



i) S^ink^fJ-k+i) = 2 iyu'^ + 16^2 - u) l^fc + ^k+i\ + 



(1) 



where N is the number of primitive cells in the 
chain; c[. ^.^ and c^.i-.a are the operators of creation 
and annihilation, respectively, of an electron with 
the spin a G {ti-H a-t the interstitial position {k,i); 
nk,i;G = c|, j^.^Ck,i-a is the operator of electron num- 
ber; fik is the z-component of the spin-1/2 operator; 
Sk,i — ('T-fc,i;t ~"fc,i;i)/2 is the z-compoucnt of the op- 
erator of total electron spin at the interstitial position 
(fc,i); t is the transfer integral; U is the magnitude 
of the on-site electron-electron attraction {U >0); h 
and I2 are the parameters of the Ising interaction for 
the bonds along the diamond sides (this interaction 
being identical only for collinear bonds) (Fig. 1); and 
he and hi are magnetic fields that act on electron spins 
and Ising spins, respectively. It should be noted that 
Hamiltonian ([T|) also corresponds to a simple Ising- 
Hubbard chain (the nodes and interstitial positions 
are aligned), in which the Ising spin ^k is coupled 
with the first, /i, and second, /2, neighbors. 



+ Aal^tfe -Mfe+i| - -^{^^k + t^k+i), 



£:5(^fc,Aifc+i) = -- (ViPTm^ + u^ |Aife + /i/c+i|+ 



Aal^fe -A^fe+il - -^if^k +A*fe+i) 



£6{pk,^J■k+l) = -u - -^{fik + fj.k+i), (2) 

where A3 < Ai < A2 are the eigenvalues of the matrix 



i^ 

^ = I ^^ 2i 
2t ~U) 



From this spectrum, the ground state and the param- 
eters of a decoration-iteration transformation can be 
determined 1241 . 
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3. Numerical Results and Their Discussion 

Consider the properties of our chain in the case of 
the antiferromagnetic Ising interaction (/^ > 0) , when 
there exists a geometrical frustration in it. The mag- 
netic fields are put to be identical, h = h\ = h^. 
Without any loss of generality, we adopt that /i > I2 
and introduce the parameter A/ — /i — I2, similar 
to what was done in work [16|. Let us pass to the 
dimensionless parameters. 



t 



It will be recalled that the parameter A/ e [0,1] char- 
acterizes the degree of asymmetry for the Ising inter- 
action along the diamond sides |24| . 

Consider firstly the properties of the system in the 
ground state. The latter corresponds to the minimum 
energy in spectrum ([2]) for all possible configurations 
pLk and /ife+i. Depending on the parameters i, [/, A/, 
and /i, the system concerned can be in four ground 
states - similarly to what take place at t/ ^ |l6lJ24 1 - 
namely, the saturated paramagnetic (SPA), ferrimag- 
netic (FRI), nonsaturated paramagnetic (UFA), and 
nodal antiferromagnetic (NAF) states. The energies 
of those states per cell are |2J| 
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Fig. 2. Phase diagram (A/, h) for the ground state. Transi- 
tion hnes for two sets of i- and (/-values are shown. The results 
obtained for U = coincide with those of work llal 



where \+)k = \ t)fcj and \—)k = \ \)k describe the 
state of spin ^k, \ t,t)fe,i;fe,2 = 4,i;t4.2:tl0), and the 
notation (— )" means the sign of (—1)". The other 
notations are 



^V- 



-At 



(it,;) + i;,t))+s^+(in,o) + |o,n)), 



^SPA — 1 7: ^, i^FRI — — 1 H — 



*^^=^5^^i t,i)+^5^^i iA)+B+-{\ n,o)+|o,n)), 

where 



i^UPA ^ -7:\U 



+ ^1/2 + mA 



h 
2' 



£1 



NAF 



= A3, 



where A^ are the eigenvalues of the matrix C = C/Ii . 
The wave functions of those states are 12^] 



N 



|SPA) = ni+)Jt,t)M;M' 

N 

iFM)=ni-)^it.t)M 



fc,2' 



fe=l 



N 



n \H'\ 



|NAF) = 



fe=i 

N 



* 



fc,2 ■ 



(-)"(-)' 



fc,l;fc,2 



n |(-)^-+^), 



,/c=l 



* 



k + l,_-.k + 2 



i-r^'i-) 



k,l;k,2 
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t,;) 



Cfc,l;t''fe,2;i 



|o), U,t) 



''fe,l;i'^fc,2;t 



|0), 



I n,0) = cl,,^cl,JO), |0,n) = 4.,2;t4,2;JO)> 



and the quantities A^ , 



-°5 I 



A^+, and B+- 



are obtained from the corresponding coefficients in 
work [2J| by replacing U ^^ —U . 

Let us consider the phase diagram for the ground 
state, {AI,h), which depends on the parameters t 
and U. One of the possible phase diagrams (A/, h) 
is depicted in Fig. 2. The FRI and UPA states in it 
are separated by the transition line A/f|u — Tit, U), 

where f{i, U) = 2 - U - ^/TpTi^. The criti- 
cal point between the FRI and NAF states in the 
zero field, A/f.n, is determined from the equation 

A/f.N — 2 + 2fNAF- 

The topology of the phase diagram (A/, h) for the 
ground state belongs to one of three types, depending 
on the parameters t and U . The topology type is 
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Fig. 3. Topological diagram (t, U) for the phase diagram of 
the ground state, {AI,h), with "equitopological" lines. Bold 
lines separate the ranges of three typical topologies designated 
by the corresponding figures 

determined by the parameter T{t, U). The first type 
of the phase diagram topology is realized at 1 < T. 
In this case, it looks as a part of Fig. 2 within the 
interval [0, A/f|u] [16|,|2J|. The second topology type 
is reahzed at < T < 1 (see Fig. 2). The third 
topology type is realized at T < Q. In this case, 
the phase diagram (A/, h) looks like a part of Fig. 2 
within the interval [A/pm, 1], but now the line of the 
NAF -H- UPA transition always begins at the point 
(0,0) MM- 

It is convenient to represent the dependence of 
the (A/, h) phase diagram topology on the param- 
eters t and U in the form of a topological diagram 
{t,U), which is shown in Fig. 3. "Equitopological" 
lines in this diagram are described by the equation 
T{t, U) =const. The changes of U and t along the 
"equitopological" line are reflected in the phase di- 
agram (A/, h) as a displacement of only those lines 
that bound the range of the ground NAF state, which 
is demonstrated in Fig. 2. In particular, as the pa- 
rameter U grows, the range of the ground NAF state 
in it decreases. 

Now consider the influence of the on-site electron- 
electron attraction on the magnetization, magnetic 
susceptibility, and heat capacity in the regime 
T(<, tj) — const, i.e. along an "equitopological" 
line. Numerical calculations of those characteristics 
were carried out for a number of points (i, U) along 
the "equitopological" line that belongs to region 2 in 
Fig. 3 and passes through the points (0.375,0) and 
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Fig. 4- Temperature dependences of the zero-field heat ca- 
pacity at various A/: (a) A7 < A/p n and (6) A/ > A/p jm. 
Solid and dashed curves correspond to the results obtained for 
t = 0.15 and U = 0.63. Dotted curves show the results ob- 
tained for i = 0.375 and U = 0, which were exhibited in work 
[iql except for the curve for A/ = 0.4 



(0.15,0.63), for which the phase diagram {AI,h) of 
the ground state is exhibited in Fig. 2. 

A comparison of the results obtained for the mag- 
netization and the magnetic susceptibility at various 
points on the "equitopological" line shows that the 
field and temperature curves of magnetization and 
the temperature curve of magnetic susceptibility in 
the zero field shift at strengthening the attraction, 
similarly to what takes place at weakening the on-site 
electron-electron repulsion |2J]. In particular, if the 
attraction becomes stronger, the temperature curves 
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of total and electron magnetizations shift downward, 
and the temperature curve of magnetization for Ising 
spins shifts upward. 

The modification of the temperature dependence 
of the zero-field heat capacity under the influence of 
the attraction is shown in Fig. 4. In the absence of 
attraction ((7 = 0), the temperature curve of the heat 
capacity has the main and low-temperature maxima 
in a wide range of A/ [16] . In a certain vicinity of the 
critical point A/f.n, it has another low-temperature 
maximum, which is the closest to the zero temper- 
ature [2J|. While the attraction grows to a defi- 
nite value of about 0.42, the main maximum shifts 
toward lower temperatures. However, if the attrac- 
tion grows further, the maximum shifts back toward 
higher temperatures. As a result, the main maxi- 
mum of the curves that correspond to small A/ can 
merge with the low-temperature maximum (Fig. 4). 
Moreover, starting from U ~ 0.5, there emerges an 
additional maximum in the temperature dependence 
of the heat capacity, which is located between the 
low- temperature and main maxima (Fig. 4). At first, 
this additional maximum exists within a small inter- 
val for A/ in the range of the ground NAF state. As 
the attraction becomes stronger, this interval broad- 
ens and surrounds the critical point A/f.n (Fig- 4). 
The modification in the temperature dependence of 
the zero-field heat capacity induced by the attraction 
growth is associated with changes in spectrum ([2]) of 
the cell Hamiltonian T-Lk- Namely, as the parameter 
U increases, the energies £q, f4(A2), and f3(Ai) in 
this spectrum decrease. 

It is also worth noting that the temperature curve 
of the heat capacity for U — 0.63 has an additional 
low-temperature maximum in a very close-much nar- 
rower than at U — 0-vicinity to the critical point 
A/p.N- This maximum, which is the nearest to the 
zero temperature, was described in detail in work |24| . 

4. Conclusions 

In this work, the properties of an asymmetric dia- 
mond Ising-Hubbard chain with the on-site electron- 
electron attraction are studied in the ground state 
at finite temperatures. It is an example of the ex- 
act solution obtained with the use of the decoration- 
iteration transformation. In the case of the antiferro- 
magnetic Ising interaction, when there is a geometri- 
cal frustration in the chain, the influence of the on- 
site electron-electron attraction on the ground state. 



the field and temperature dependences of the mag- 
netization, and the temperature dependences of the 
zero-field magnetic susceptibility and the heat capac- 
ity are studied. 

The phase diagram for the ground state is plot- 
ted in the plane {AI,h). A modifications of this 
phase diagram induced by the transfer integral and 
the on-site attraction is represented in the form of a 
topological diagram (i, U). It is shown that strength- 
ening the attraction along the "equitopological" line 
T{t, U) ~ const is reflected in the phase diagram 
(A/, h) as a displacement of the ground NAF state 
boundaries so that the corresponding confined area 
becomes smaller. 

As the on-site electron-electron attraction becomes 
stronger along the "equitopological" line, the temper- 
ature curves of the magnetization, at various fields, 
and the zero- field magnetic susceptibility shift, as it 
was in the case of weakening the on-site electron- 
electron repulsion [24!| . The temperature dependence 
of the zero-field heat capacity in a certain range 
(i, U , A/) has an additional maximum between the 
main and low-temperature maxima. 

The results of this work also correspond to a simple 
Ising-Hubbard chain, in which the nodes and the in- 
terstitial positions are aligned, and the nodal Ising 
spin is coupled with those of the first and second 
neighbors. 
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Translated from Ukrainian by O.I. Voitenko 

B.M. JlicHuu 

ACHMETPHHHHH POMBIHHHH JIAHll,IO>KOK 
ISIHrA-rABBAPAA 3 nPHT5irAHH51M 

P e 3 K) M e 

Po3rji5iHyTO ocHOBHiiH CTaH i TepMOfliiHaMinni BjiacTiiBOCTi 
aciiMeTpHHHoro poM6i'^Horo jiaHLi,iO>KKa Isinra— ra66apfla 3 

O.Z],HOn,eHTpOBHM ejieKTpOH-ejieKTpOHHHM npiITHraHHHM, HKHli 

e TOMHO po3B'5i3yBaHHM 3a flonoMoroio fleKopai^iHHO-iTepau;iH- 
Horo nepeTBopeHHH. Y BHna^Ky aHTiicbepoMarHiTHOi B3aeMO- 
flii' Isinra BHBMeno brjimb Li;boro npHTHraHHH na ochobhhh 
CTaH i TeMnepaTypny 3ajTe:>KHicTb HaMarHineHOCTi, MarniTHOi 

CnpHilHHTJlHBOCTi Ta TenjlOCMHOCTi. 
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